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Abstract. - Universality of strongly interacting fermions is a topic of great interest in diverse 
fields. Here we investigate theoretically the universal dynamic density response of resonantly 
interacting ultracold Fermi atoms in the limit of either high temperature or large frequency: (1) At 
high temperature, we use quantum virial expansion to derive universal, non-perturbative expansion 
functions of dynamic structure factor; (2) At large momentum, we identify that the second-order 
expansion function gives the Wilson coefficient used in the operator expansion product method. 
The dynamic structure factor is therefore determined by its second-order expansion function with 
an overall normalization factor given by Tan's contact parameter. We show that the spin parallel 
and antiparallel dynamic structure factors have respectively a tail of the form ~ ±tj~^''^ for 



CO, decaying slower than the total dynamic structure factor found previously 
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Our 



predictions for dynamic structure factor at high temperature or large frequency are testable using 
Bragg spectroscopy for ultracold atomic Fermi gases. 
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' Introduction. — The study of strongly interacting 
fermions has brought together very different areas of 
physics - neutron stars, quark-gluon plasmas, high tem- 
perature superconductors, and cold atoms - which, at first 
^ight, have little in common. There is, however, an im- 
portant generic idea of fermionic universality behind [I]: 
ill strongly interacting, dilute Fermi gases should behave 
identically, depending only on scaling factors equal to the 
dverage particle separation and/or thermal wavelength, 
but not on the details of the interaction. Recent manipu- 
lation of ultracold Fermi gases of ^Li and atoms near 
a broad collisional (Feshbach) resonance provides an ideal 
avenue to understand this fermionic universality PHS] . To 
date, universal thermodynamics of strongly interacting 
fermions has been clearly demonstrated 0, by measuring 
the static equation of state. The purpose of this Letter 
is to show that universality appears in dynamical proper- 
ties as well. We derive exact results for universal dynamic 
structure factor (DSF) at high temperatures or at large 
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momenta and frequencies. 

Exact results are very valuable for strongly interacting 
fermions, due to their non-perturbative, strongly corre- 
lated nature. There is no small interaction parameter to 
control the accuracy of theories. In specific cases, ab-initio 
calculations are possible using Monte Carlo methods [THH] ■ 
However, in general this approach suffers from the fermion 
sign problem [TU]. In 2008, Tan derived a set of exact, uni- 
versal relations for two-component (spin-1/2) Fermi gases 
with a large s-wave scattering length a These uni- 

versal relations all involve a many-body parameter called 
the contact T, which measures the density of pairs within 
short distances. Tan's relations can be understood us- 
ing the short-distance and/or short-time operator produc- 
tion expansion (OPE) method [12lI14j . which separates 
in a natural way few-body from many-body physics. For 
the "interaction" DSF, l\S„„' (q,w,T) = S*^^/ (q,w,T) - 
5*^^^ (q, Lo,T), the OPE predicts that (q -H' oo and w — > oo) 
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where {cr,a') =t,i denote the spin, S^cr' and S'^^l are 
respectively the DSF of interacting and non-interacting 
Fermi gases at the same chemical potential /i and temper- 
ature T, and Waa' are the temperature-dependent Wil- 
son coefficients that rely only on few-body physics. At 
high temperatures, quantum virial expansion provides an- 
other rigorous means to bridge few-body and many-body 
physics |15H18| . It was shown that static thermodynamic 
properties of a strongly correlated Fermi gas can be ex- 
panded non-perturbatively in fugacity using some univer- 
sal, temperature- independent virial coefficients jSl lTTllTS] . 
which are exactly calculable from few-fermion solutions. 
Both OPE and virial expansion give useful insight into 
the challenging many-body problem. However, their con- 
nection is yet to be understood. 

In this Letter, we investigate theoretically the univer- 
sal dynamic properties of a strongly correlated Fermi gas 
in the limit of either high temperature or large momen- 
tum/frequency. In the former limit, we show that the 
dynamic structure factor can be virial expanded in fu- 
gactiy, using some universal, temperature-implicit virial 
expansion functions. We derive, for the first time, these 
universal virial expansion functions for spin parallel and 
antiparallel DSFs of a homogeneous Fermi gas in the res- 
onance (unitarity) limit, where the scattering length di- 
verges (a — ±00). In the latter limit of large momentum, 
we show that the Wilson coefficient in Eq. ([T|) is given by 
the second-order expansion function. Therefore, the large 
momentum DSF is universally determined by the second 
order virial expansion, together with a many-body pref- 
actor - the contact. Our results can be easily examined 
using Bragg spectroscopy for ultracold Fermi gases of ^Li 
or atoms [H]. 

Universal expansion function of DSF. Virial 
expansion is a powerful tool for studying the high- 
temperature properties of ultracold atomic Fermi gases 
[TT] . It expresses any physical quantities of interest as an 
expansion in fugacity with some coefficients or functions, 
since the fugacity z = exp{fj./kBT) ^ 1 is a controllable 
small parameter at high temperatures. For the interaction 
DSF, the expansion is given by fIU\ 



(2) 



where AiS'crcr',)i(q, w, T) {n = 2,3,---) is the n-th expan- 
sion function. The expansion is non-perturbative as the 
few-body problem could be solved exactly, no matter how 
strong the interaction. Virial expansion is anticipated to 
work for temperatures down to the superfluid transition, 
although a nontrivial resummation of expansion terms 
may be required if z 1. In Ref. [20], the lowest second 
order expansion function AS'^^"^'^ for a trapped Fermi gas 
was calculated using two-fermion solutions in traps. How- 
ever, the universal aspect of expansion functions was not 
realized. As a result, for different temperatures/momenta 
the expansion functions had to be re-calculated. 
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Fig. 1: (Color on-line) Universal second order expansion func- 
tion of DSF ai q = 1/3, 1, and 3. The inset shows the rapid 
convergence of A5'2(g, (i) at small fiujo / {ksT) (thick lines) and, 
AS'ft,2 (thin solid line) and A5-|-4,,2 (thin dashed line) ai q = 1. 



Here we consider the expansion functions of a homo- 
geneous Fermi gas in the unitarity limit and emphasize 
their universal aspect, which is not known so far. As the 
scattering length diverges, all microscopic scales of the in- 
teraction are absent [1]. For this few-body problem, the 
only energy scale is /c^T and length scale is the thermal 
wavelength A = /i/(27rmfcBT)^/^. Dimensional analysis 
leads to, 



A5,,.,„(q,w,T) 



V 



AS'acr',n(<7,w), 



(3) 



where V is the volume, q = [?i^q^/(2mA;BT)]^/^, Cj = 
fujj / {kBT), and AS'o-o-'.n is a dimensionless expansion func- 
tion. The temperature T is now implicit in the variables 
q and Cj. This universal form implies a simple relation 
between the trapped and homogeneous expansion func- 
tion. In a shallow harmonic trap, VTrap(r) = TO(a;^a;^ -I- 
ujyV'^ +w1z'^)/2, where cjq = {uJxUJyUJzY^^ 0, the sys- 
tem may be viewed as a collection of many cells with a 
local chemical potential /x(r) = ^ — ^rap(r) and fugacity 
z{r) — zexp[— VTrap(r)/fcBr], so that the trapped DSF is 

given by AsZ".""] (q,a;,r) = ^ dv[AS,,, (q, T, r) /F]. 
Owing to the universal q- and w-dependence in the ex- 
pansion functions, the spatial integration can be easily 
performed, giving rise to 



ASaa' ,n{q,C]) 



,3/2 



^A5(^rj)(q,c.,T). (4) 



The (non-universal) correction to the above local density 
approximation is at the order of 0[(?ia;o)^/(fc_BT)^]. Eq. 
(|4]) is vitally important because the calculation of expan- 
sion functions in harmonic traps is much easier than in 
free space. 

Fig. 1 reports the homogeneous expansion function 
A52 = 2[AS'-|-|-,2 + AS^i^2] at three different momenta. 
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using AS'^"^r2^^ ^^f- EO] as the input. One observes 
a quasielastic peak at w = 9^/2 or a; = ?iq^/{4m), as a 
result of the formation of fermionic pairs. We note that 
the third expansion function A^o-<t'.3 or AS'^^r'^^^ can also 
be calculated straightforwardly using exact three-fcrmion 
solutions P^ . 

We may derive sum rules that constrain the universal 
expansion functions, using the well-known /-sum rules sat- 
isfied by DSF: toS^ici, w, T)duj = Nq^/{4m) ^ and 
J^^ ojSfi{ci,uj,T)duj = ED. The latter immediately 
leads to 

f +00 

a;A5i-4„„(g, = 0. 



(5) 



On the other hand, virial expansion of the total number 
of fermions N implies that 



+ CXJ 



J)AS'^^^„(g,a))cfa; = nq^Abn, 



(6) 



where A6„ is the n-th virial coefficient and in the unitarity 
limit A62 = 1/^2 and A63 ~ -0.355 [5j[l7]. 

At large momentum, the spin-antiparallel static 
structure factor satisfies the Tan relation j21| . 
/ 5-t-4.(q, w, T)du! ~ I/{8hq). This indicates a virial expan- 
sion of the contact: I = l6n^V [z^C2 + z^cz + • • •] /A"', 
where the contact coefficients c„ are given by, 

A%.„(q»l)= / A%,„(g,w)dw = — ^. (7) 

J —oc 9 

The expansion of the contact was alternatively obtained 
using Tan's adiabatic sweep relation |18j . In the unitarity 
limit, it was shown that C2 = I/tt and C3 ~ —0.141 [TH]. In 
the same limit of large momentum, the spin-parallel static 
structure factor is nearly unity so that / 5'tt(q, ui,T)dui ~ 
N/{2K) [inilll]- This leads to 



A%,„(g> 1) 



+ 00 



(8) 



For the second expansion function, AS'o-cr',2, we have 
checked numerically that all the above mentioned sum 
rules arc strictly satisfied. 

Wilson coefficient of DSF. — We now turn to the 
large momentum/frequency limits, where the OPE Eq. (1) 
is assumed to be applicable. It is clear from the equation 
that the Wilson coefficient determines entirely the DSF at 
large (q, oj) as far as the many-body contact is known. 

At T~0, the Wilson coefficient W^a' can be calcu- 
lated using Feynman diagrams [13] for dynamic suscepti- 
bility Xt7c7'(r, t) = - (T^/5^(r,r)/5^'(O,0)), as S'<^<^'(q,w) 
-lmXaa'{(\,uj)/[K{l-e-^'^'^'''^)]. In the limit of (q,w) ^ 
00, the contributing diagrams to Xcrcr'(q, ^) are sketched 
in Fig. 2 [T3]. Diagrammatically, the contact may be 
identified as the vertex function at short distance and 
time [HISS]: I = -m^T{r = 0,t = 0-)/?i'^. Therefore, 
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Fig. 2: (Color on-line) Diagrammatic contributions to the 
interaction dynamic susceptibility. The self-energy (S) and 
Maki- Thompson (MT) diagrams contribute to Ax^^{r,T) and 
Ax^l{r,T), respectively, while the Aslamazo-Larkin diagram 
(AL) contributes to both. The shadow in the diagrams repre- 
sents the contact I. The crossed part in the diagram (AL) is 
the vertex. 



in the diagrams the shadow part of the vertex function 
r(r = 0,r = 0~) represents the contact I. These dia- 
grams are well-known in condensed matter community. 
In the context of calculating the change in conductiv- 
ity due to conducting fluctuations, the last two diagrams 
are called the Maki-Thompson (MT) [55] and Aslamazov- 
Larkin (AL) contributions [U] respectively, while first di- 
agram gives the self-energy (S) correction. At zero tem- 
perature, we calculate these diagrams in vacuum at ^ = 

ifs - /AL)/(47r2V^c^3/2) and 



and obtain that 
Wfr° = (f^'T - /AL)/(47r2V^c^3/2), where. 
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X = q-^ / {2mhuj) , and 8 is the step function. These 
results agree with the previous calculations by Son and 
Thompson [14], although there the spin parallel and an- 
tiparallel DSFs were not treated separately. At small 
q^/w, we find that the spin parallel and antiparallel DSFs 
have the tail 



127r2m3/2w5/2 ■ 



(9) 



This prediction shows that for w — >■ 00 the spin dependent 
DSFs decay an order slower in magnitude than the total 
DSF. The latter is proportional to 5^/0;^/^, as shown in 
Refs. [2] and [55]. The faster decay in the total dynamic 
structure factor is due to the cancellation of the leading 
terms in W^=° and W^^=°. 

It is not clear how to obtain the finite temperature Wil- 
son coefficient using diagrammatic technique, since the fi- 
nite temperature diagrams involve many-body process in 
medium. However, Eq. (1) gives the hint. It has a strong 
constraint on the expansion functions of DSF. As W^a' 
involves only the few-body physics and hence does not 
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Fig. 3: (Color on-line) f„„, = ^/m?i^J^''^(2^/^2) AS'^^/_2 at g = 
3, 5, and 10. With increasing momentum and/or frequency, 
/ctct' approaches smoothly to the T = Q result by Son and 
Thompson [Tl] . 



contain the fugacity z, a count of tlie term z^' 
sides of Eq. (1) leads to 



both 



(c„/c2)A5,,,,2(q,c^,T), (10) 
iz^ /I2) I^^S^a' ,2 (q,c^,T),(ll) 



where I2 = 2;^167r^yc2/A^ is the contact up to the second 
order expansion [18]. Therefore, the Wilson coefficient is 
given by the second expansion function. This result is ob- 
tained by applying the OPE and virial expansion method. 
As a result, in principle it should be valid at temperatures 
above the superfluid transition. However, we may expect 
that it holds at all temperatures, as both the Wilson coef- 
ficient and second expansion function are irrelevant to the 
many-body pairing in the supcrfiuid phase. The many- 
body effect enters through the many-body parameter of 
contact only. 

Eq. pT|) is the main result of this Letter. At a first 
glance, it may be a suprising result. However, it could be 
understood from the proportionality shown in Eq. (jlO[) . 
On the other hand, a rigorous proof of Eq. pU)) at large 
q and w justifies the use of the OPE method. We note 
that, for the two-component Fermi gas, our result shows 
that the Wilson coefficient is determined by the two-body 
physics solely. This is in agreement with previous work. 
For example, in Ref. [13] a two-body scattering ampli- 
tude was used to calculate the zero temperature Wilson 
coefficient for the rf-spectroscopy of a strongly interact- 
ing Fermi gas. However, in general cases where three or 
four-body physics come into play, we anticipate that the 
Wilson coefficient should be related to the higher order 
virial expansion function. In that case, new universal re- 
lations with new many-body "contact" parameters would 
appear. 




Fig. 4: (Color on-line) Universal total DSF (a) and spin- 
antiparallel (b) DSF at T /Tf = 0.225, 0.5, 1, and 2, calculated 
using Eq. (|12|) . We use a pair fluctuation theory to determine 
the fugacity z and contact T [18j . 



In Fig. 3 we check the vafidity of Eq. ([TT]) at T 0, 
by calculating y/nihuj'^/'^{z'^ /l2)lS.Saa' ,2 at different mo- 
menta. With decreasing temperature T or increasing 
q oc q/Vr, it approaches gradually to y/mTwj^/'^Wj^'^ . 
This confirms numerically that Eq. pT]) holds at zero 
temperature. Moreover, at high temperatures where the 
fugacity is small, to a good approximation we have the in- 
teraction DSF ASaa' — z^ASaa' ,2- As the contact I ~ Z2 
at high T, it is trivial to confirm Eq. ([TT]). Note that, in 
the limit of large frequency, the tail w"^/^ of and 
W^.j^^'^ is fairly evident in the second order virial expansion 
functions. 

Universal DSF at large (q, w). — In this limit, using 
Eqs. ([!]) and pT]) the DSF is approximated by 



This approximate DSF should be quantitatively accurate 
for sufficiently large momentum and frequency. It holds at 
all temperatures and the many-body effect is respected by 
the contact. However, the momentum q should be larger 
than a critical momentum q^ ^ max{X~^ ,kp) in order 
to validate the use of the OPE equation ([TJ. Here, kp 
is the Fermi wavevector. Quantitatively, an estimate of qc 
requires the calculation of AS^a' ,n>2 and the examination 
of Eq. ([TU| . We note that Eq. may not be reliable at 
small frequency, w ~ 0. As a result, the structure factor 
sum-rules may not be strictly satisfied. We note also that, 
in the limit of high temperatures where the contact X ~ T2, 
the prefactor of the ASaa',2 term is z^. Therefore, the 
approximate DSF reduces back to the virial expansion up 
to the second order. In this high-T limit, Eq. ([12]) is valid 
for arbitrary q and lu. 

We present in Fig. 4 the temperature dependence of 
the approximation DSF of a normal, homogeneous uni- 
tary Fermi gas at q = bkp, by assuming that qc ^ 5kp. 
The many-body contact and fugacity are calculated by us- 
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ing a strong-coupling pair fluctuation theory [TB], which 
is shown to be accurate for describing the unitary equa- 
tion of state. Close to the superfluid transition temper- 
ature, a quasielastic peak clearly emerges at the recoil 
energy for pairs, uj = ?iq^/(4m), in agreement with the 
low-temperature experimental observation |19| . 

Observation of universal DSF at large frequency. 

— Eq. ([T^ indicates that at large (q, w), the inter- 
action DSF of a unitary Fermi gas depends on the re- 
duced moment q = qX/^/in and reduced frequency uj 
only. This universal dependence could be examined using 
large- momentum Bragg spectroscopy |191I21| . with vary- 
ing momentum and temperature while keeping q invariant. 
One can also extract experimentally the universal second 
expansion function AS'o-cr'.2(<7, since the contact I can 
be determined independently using the /-sum rule |21j . 
These predictions break down below the critical momen- 
tum qc- Note that, by tuning the transferred momentum 
in Bragg beams, the value of qc 3> raax{\^^ ,kp) might 
be determined experimentally. 

Conclusion. — We have studied the finite tempera- 
ture dynamic structure factor of a homogeneous unitary 
Fermi gas, using quantum virial expansion and operator 
product expansion. The universal second order expansion 
function has been calculated and related to the Wilson co- 
efficient at large momentum q. As a result, in that limit 
the thermal wavelength A becomes the only length scale 
and the interaction dynamic structure factor should de- 
pend universally on a reduced momentum gA/\/47r- We 
have proposed that Bragg spectroscopy with large trans- 
ferred momentum should be able to confirm this universal 
dependence. Our results can be extended to other dynam- 
ical properties of a strongly correlated Fermi gas, such as 
the rf-spcctrum and single-particle spectral function. 

It is an interesting chanllenge to derive from three- and 
four-fermion solutions j221l29] new universal relations in- 
volving a many-body parameter like Tan's contact. The 
determination of Wilson coefficients in that case should 
be difficult. Our method of calculating higher-order virial 
expansion function would give the most natural and con- 
venient way to obtain the Wilson coefficient. 
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